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Abstract Our previous semi-analytic treatment of /- and p-mode multiple 
scattering from ensembles of thin flux tubes (Hanson and Gaily, Astrophys. J. 
781, 125; 791, 129, 2014) is extended by allowing both sausage and kink waves 
to freely escape at the top of the model using a radiative boundary condition 
there. As expected, this additional avenue of escape, supplementing downward 
loss into the deep solar interior, results in substantially greater absorption of 
incident /- and p-modes. However, less intuitively, it also yields mildly to sub¬ 
stantially smaller phase shifts in waves emerging from the ensemble. This may 
have implications for the interpretation of seismic data for solar plage regions, 
and in particular their small measured phase shifts. 
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1. Introduction 


The interaction of the Sun’s /- and p-mode oscillations with large collections 
of thin magnetic-flux tubes that make up plage and network is important for 
our understanding of the broader seismology of solar activity. Individual tubes 
absorb, phase shift, and scatter incident waves, but close assemblies of tubes 
provide a multiple-scattering regime and collective behaviour that makes them 
more than just a summation of their parts. 

The historical context is addressed in detail by Hanson and Gaily ( 2014al f9l. 
Briefly, early modelling of multiple scattering in ensembles of flux tubes ignored 
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gravitational stratification due t o the considerable mathema tical c omplications 
that aris e (Bo gdan and Zweibel, IQSSt Zw eibel and Bogdan, EHi Bogdan and 
Zweibel, Ryutova and Priest, 

Techniques for dealing with scattering from single tub es in s tratified atmo¬ 
spheres w ere de veloped by, for ex ample, Hanasoge et al. (1200811 a nd Hiii dman 
and Jain ( 2012h . Jain et al. ( 2009), Jai n, Gascoyne, and Hindman ( 2011a ). and 
Jain, Gascoyne, and Hindman ( 201Hj) considered multiple tubes, but did not 
allow for multiple scattering, thereby treating the ensembles as simply the sum 
of the individual tub es. Th e true multiple-scattering regime was attacked by 
Hanasoge and Cally ( 2009h for t he cas e of two adjacent tubes, adapting the 
formalism of Kagemoto and Yue ( 19861) from the engineering literature. Only 
kink waves were modelled, as it was not then clear how to implement the proper 
side boundary cond itions for sa usage modes. The role of the acoustic jacket 
(Bogdan and Callv. [l995l : Cally. [2013 ) was emphasised in mediating near-held 
coupling between tubes. 

The proper approach to sid e bou ndary conditions for saus age m odes was 
resolved by Andries and Cally (120111) a nd Hindman and Jain (|2012l) . and im- 
plemen ted by Hanson and Cally ( 2014al) for two tubes, and Hanson and Cally 
( 2014bl) for multiple tubes, again using the Kagemoto and Yue ( 1986h formalism. 
A slight modihcation was made to correct an unp hysica l symmetry breaking that 
appeared in the results of Hanasoge and Cally (l2009l) . Direct numerical simu¬ 
lations (not using the thin-tube appro ximati on) yielded gener al conhrmation of 
the results (Felipe, Crouch, and Birch. [201^ Daiffallah. 12013 ). 

Wave absor ption and phase advance ment h a ve been measured in solar plage 
(Braun, Il995l : Chou, Sun, and Chang, l2000l) . Braun reported no meas urable 
phase change in waves emerging from plage, but Chou. Sun, and Chand found 
a weak positive shift. On the oth er han d, phase shifts calculated for thin-tube 
ensembles by Hanson and Cally ( 2014bl) are always negative. The likely cause 
of this discrepancy is that the ensemble models postulate field-free quiet-Sun 
plasma between the flux tubes, whereas real plage is both thermally different 
from quiet Sun and probably has a continuum of weak field in addition to the 
many “discrete” flux tubes. Clearly, continuum field speeds up transiting waves 
(which are fast MHD waves in magnetic environments) thereby advancing their 
phase, whereas multiply scattering discrete ensembles slow them down in the 
manner of a ball in a pin-ball machine, resulting in negative phase shifts. The 
competition between these two effects may in part be responsible for the small 
measured shifts from plage. 

To date, it has been customary to impose a stress-free boundary condition at 
the top of each tube, thereby precluding any upward wave losses. The purpose 
of this article is to investigate the implications of relaxing that restraint. 


2. Upwardly Radiative Boundary Conditions 

The mathematical formalism developed in Hanson and Cally ( 2014allbl ) addresses 
the absorption and scattering of waves by single and multiple flux tubes. The 
“absorbed” waves take the form of sausage or kink tube waves, with their 
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attendant acoustic jackets, that propagate freely down into the solar interior 
and are therefore lost to the observable surface oscillations. This is the nature 
of absorption in those models. However, in reality, waves may also propagate 
upward into the solar atmosphere. This was expressly prohibited in our previous 
work by imposing perfectly reflective stress- free boundary conditions at the top 
of each tube, following Bogdan et al. (1199611 . 

Th ere have been several attempts to relax this restriction. Crouch and Cally 
I I 999 II allowed the thin tube to continue upward through a crude model chro¬ 
mosphere and corona without expanding in radius and matching onto outgoing 
magneto-acoustic waves in the corona. This of course is also very unrealistic, 
since the tube will inevitably expand well beyon d the thin-tube regime very 
quickly. On the other hand, Hindman and Jain ( 200811 sidestep the issue of 
a realistic boundary condition by introducing the concept of a maximal-flux 
condition. This is applied in detail to sausage modes in single flux tubes by 
Gascoyne, Jain, and Hindman ( 2014ll . 

Ad opting the truncated adiabatic polytrope external model of Bogdan et al. 

1 1996ll . with polytropic index mp, truncation depth zq, dimensionless depth 
variable s = —zIzq (1 ^ s < 00), and assuming the thin flux tubes are of 
uniform plasma j3 (the ratio of gas to magnetic pressure), the thin tube oscillation 
equation for the sausage mode (azimuthal order m = 0) is 


w^(2mp -I- /3(mp -I- 1)) -I- 


2gs 

Zq ds"^ 


gjiTip -I- 1) 9 
ds 


Zo 


^11 


= -uj'^{mp + l){/3 + 1) 


ds 


( 1 ) 


and that of the kink mode (m = ±1) is 


Lu'^Zq 


2gs 9^ 

(1 -I- 2l3)(mp + 1) ds"^ 


9 d 1 _ 2{l + j3) 2 ^^inc 

l-b2/3asj “ l-b2/3 dx ■ 


( 2 ) 


Here and are respectively the longitudinal and transverse plasma displace¬ 
ments, g is the gravitational acceleration, and w is the circular frequency. The 
inhomogeneous term dtinc represents the external driving by /- and p-modes. 
Polytropic index nip = 3/2, corresponding to adiabatic index 7 = 5/3, is used 
throughout. 

These driven tube equations may be solved subject to appropriate boundary 
conditions at the top (s = 1) and bottom (s = 00 ) by constructing Green’s 
functions. As usual, the Green’s function is built from the homogeneous solutions 
satisfying suitable homogeneous boundary condi tions. These may be expressed 
in terms of Hankel functions (Grouch and Gallv. [l99^ . 


^/’<t(s) = s {2iy ,/^), 

</,(s) = s-^/2 fijG) (2,,^) + {2,.y/^) 


( 3 ) 

( 4 ) 
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Figure 1. The absorption and phase shift for 3 mHz /-modes {n = 0) of various cylindrical 
orders m incident upon two tubes of /3 = 1 separated by 0.25A (1.225 Mm) positioned along 
the a^-axis. The two cases studied here are the stress-free (dashed red) and the radiative (full 
black) boundaries. 


where fi = {nip — l)/2, = nipUj'^Zo/g, a represents either || (sausage) or _L 

(kink), and we have adopted the convenient parameters 


e|| 


2mp +/3(TOp + 1) (wp + 1)(1 + 2/3) 

--- and e_L =- - - 


2m. 


2mp 


( 5 ) 


The solution 'ipu is constructed to satisfy a pure radiation condition at s = c», 
based on the well-known large-argument asymptotic behaviour of the Hankel 
function. The parameter may be chosen to select an appropriate boundary 
condition at the top [s = 1]. For example, a particular choice (that is n ot use d 
here) recovers the stress-free solution in Equation (4.4) of Bogdan et al. 1 1996h . 
However, here we choose simply Ao- = 0, thereby imposing an outgoi ng rad iation 
condition at s = 1 (see also Equation (4.7) of Hindman and Jain, 2008ll . The 
identification of the ffp Hankel function as “purely upgoing” is not trivial, as 
the exponential asym ptotic behaviour does not (generally) apply there. Nev¬ 
ertheless, Cally (120121) has shown that the Hankel function at small argument 
does actually represent a unidirectional wave coupled with an unavoidable wake 

(unless n is half an odd integer) that decays algebraically in time. On that basis, 
f 2) 

we identify the Hp ’ solution as the radiation solution at the top. 

The rest of the analysis is exactly as set out by Hanson and Cally ( 2014al lbh. 


3. Results 

For conciseness, we will mainly consider incident /-modes, as these interact most 
strongly. Corresponding results for the pi mode will be presented very briefly. 

First, we address the case of two /3 =1 flux tubes (one at coordinate centre) 
separated by a quarter of an incident /-mode wavelength (0.25A). The absorption 
and phase shift for the two top boundary conditions, stress free and radiation, 
are compared in Figure [T] Absorption amn and phase shift (pmn of an incident 
wave of cylindrical degree m and radial order n are defined by Hanson and 
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Figure 2. Top: the absorption and phase shift (degrees) of 3 mHz /-modes (n = 0) of 
cylindrical order m that are incident upon an ensemble of seven i dentical tubes (/3 = 1) 
randomly situated around the origin (as in Hanson and Callv. [2014bfi . The dashed red curve 
shows the scattering properties of the stress-free boundary condition and the full black the 
transmissive boundary condition. The bottom two frames show the change in absorption and 
phase from the case of the non-interacting ensemble (multiple scattering — single scattering). 
The random positioning of the tubes breaks m —m symmetry, as is most apparent in the 
bottom-left panel. 


Cally ( 2ni4bl . Equations (4) and (5)). The /-mode corresponds to n = 0. Not 
surprisingly, the ability of the tubes to lose energy upward roughly doubles the 
absorption coefficients for both the sausage (m = 0) and kink (m = 1) modes. 
Higher-order modes (|m| > 1) are weakly seen by the system as the off-axis tube 
reacts to m = 0 and 1 components in the incident wave as experienced from 
its own coordinate centre. On the other hand, the phase shift is substantially 
reduced by the open boundary condition, with only the sausage mode retaining 
any significant shift. 

Figure [5] displays similar results for the same randomly positioned se ven-tub e 
ensemble of /3 = 1 tubes as described in Table 2 of Hanson and Cally ( 2014bli . 
Once again there is a clear difference between the stress-free and radiative top 
boundary cases, with absorption of the sausage mode more than doubled and 
the negative phase difference slightly reduced in magnitude. The lower panels 
illustrate the importance of the multiple scattering regime. Multiple scattering 
has a far larger effect for kink modes (|m| = 1) in particular where the radiative 
top boundary is applied. 

Figure 13] again extends some results from Hanson and Cally ( 2014bl) . Absorp¬ 
tion and phase shift are plotted for the sausage and kink modes, with either the 
stress-free or radiative boundary condition, for between 2 and 20 identical ran¬ 
domly placed /3 = 1 tubes as they are added, one after the other, within a radius 
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Figure 3. The absorption (left) and phase shift (degrees, right) of /-modes incident on a 
randomly distributed ensemble of tubes as a function of increasing tube number. The tubes 
are identical (^ = 1) and can scatter up to |m| = 4. The stress-fre e boun dary cases are shown 
by the black dots and red stars (identical to Hanson and Cally, l2Q14bn for the m = 1 and 
m = 0 scatter, respectively. The radiative tubes are shown by the green squares and blue 
circled crosses for the m = 1 and m = 0 waves, respectively. 
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Figure 4. The absorption (left) and phase shift (degrees, right) of pi-modes incident on a 
randomly distributed ensemble of tubes as a function of increasing tube number. The tubes 
are identical (^ = 1) and can scatter up to |m| = 4. The stress-fre e boun dary cases are shown 
by the black dots and red stars (identical to Hanson and Cally, l2Q14bn for the m = 1 and 
m = 0 scatter, respectively. The radiative tubes are shown by the green squares and blue 
circled crosses for the m = 1 and m = 0 waves, respectively. 


of 2.5 Mm from the origin. The precise details of the results depend of course on 
the random placement, but broadly absorption is significantly enhanced in all 
cases by opening up the top boundary, and increases roughly linearly with tube 
number. The phase shift is reduced in magnitude strongly for the kink wave and 
slightly for the sausage wave. 

Figure 0] repeats these calculations, but for incident pi- instead of /-modes. 
For each boundary condition, both the absorption and phase shift are reduced 
by roughly an order of magnitude compared to the /-mode case. Once again, net 
absorption is increased by opening the flux-tube tops to wave egress, especially 
for the kink mode. The phase shifts roughly halve in both cases. 

Finally, in Figure [S] we again compare the 2-20 randomly placed tubes, but 
for 13 = 1, 0.1, and 10. Absorption of sausage modes (top-left panel) is most 
pronounced for the strongest tubes {13 = 0.1), though note that through hap- 
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Figure 5. The absorption (left) and phase shift (degrees, right) of the m = 0 (top) and m = 1 
(bottom) incident and outgoing /-mode waves interacting with a random ensemble of growing 
tube number of identical properties (/? = 1: filled black circles; p = 0.1: filled red triangles; 
/3 = 10: open green squares) with the radiative (points with lines) and stress free (points only) 
boundary conditions. 


penstance, tubes 17 and 18 produce a particularly strong increase in absorption 
for /3 = 0.1. The phase shift is broadly reduced by the radiative boundary, 
especially for the weak tubes (3 = 10. 

For the kink modes (bottom panels of Figure [S]), the weak-tube radiative case 
is clearly dominant in absorption, and most muted in phase shift. This compares 
with the stress-free boundary case, where /3 = 10 has the greatest phase shift. 


4. Conclusions 

Broadly, allowing tube waves to escape at the top into a notional solar atmo¬ 
sphere increases the total absorption of an ensemble of thin flux tubes, as we 
might expect. Details vary, depending on plasma /3, tube arrangement, and mode 
(sausage or kink, / or pi), but the trend is consistent. The most spectacular 
increase is seen for both m = 0 and m = 1 /-modes incident on a collection 
of up to 20 weak tubes (/3 = 10), with strong tubes’ absorption increased only 
marginally. 

The more surprising result though is that the wave losses at the top of the 
tubes generally reduce the overall phase shift produced by the ensemble, partic¬ 
ularly for weak tubes. This applies to both /- and pi-modes. We do not have a 
simple physical explanation for this, but it is significant for the interpretation 
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of the seismology of plage, and potentially provides a means of estimating the 
strengths of constituent flux tubes on average. 

The greatly reduced absorption of p -modes relative to /-modes was docu¬ 
mented in Figure 5 of Hanson and Cally (l2014bll for the stress-free top boundary, 
and carries over to the radiative case. Nevertheless, the seemingly linear increase 
in absorption with number of tubes suggests that even p-modes could measurably 
suffer absorption in plage, especially if they have lossy tops. This reduced effect 
with increasing radial order is strong evidence that the scattering and absorption 
are closely tied to the surface. 

Of course, in reality the behaviour of tube waves when they reach the photo¬ 
sphere will not be as simple as suggested by the radiative boundary condition. It 
will depend on the rapid flux tube expansion expected there, field connectivity, 
frequency, and many other details of the physics. Our purpose here has merely 
been to explore the consequences of the most extreme behaviour, total loss to 
the atmosphere of tube waves reaching the surface. This gives some insight into 
the range of behaviours that we might expect to see in the real Sun. 
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